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A WING  WITH  A SMALL  ASPECT  RATIO  IN  A RESTRICTED  FLOW  OF  A NON- 
VISCOUS  FLUID 

V.  I.  Kholyavko,  Yu.  F.  Usik 

The  aerodynamic  characteristics  of  a thin  flat  wing  which  is 
moving  close  to  a solid  or  free  surface  (hydrofoil)  are  determined. 
Here,  we  use  the  theory  of  a thin  body  in  accordance  with  which 
the  three-dimensional  flow  around  a thin  body  which  is  elongated 
in  the  direction  of  movement  is  replaced  by  a two-dimensional  flow 
in  transverse  planes.  In  the  case  under  consideration  the  problem 
is  reduced  to  a study  of  the  movement  of  a flat  plate  close  to  a 
surface  of  separation  (Fig.  1). 

Let  a flat  wing  of  small  aspect  ratio  the  maximum  span  of. 
which  coincides  with  the  trailing  edge  move  at  a small  angle  of 
attack  a with  constant  speed  V^.  During  time  dt  a section  of  the 
wing  dx.^  = V^dt  passes  through  a fixed  transversed  plane  ZOY . In 
an  equivalent  two-dimensional  problem  this  movement  of  the  wing 
corresponds  to  the  vertical  displacement  of  a flat  plate  which  is 
moving  with  constant  velocity  VQ  = V^a  by  the  value  VQdt  and  to 
a change  in  the  width  of  the  plate  from  Ztx.^)  to 

'<*.)  + - /(*,)  + 

where  - the  local  wing  semispan. 


The  flow  in  plane  ZOY  which  is  caused  by  the  vertical  displace- 
ment of  the  plate  and  the  change  in  its  width  leads  to  a change  in 
the  connected  mass  of  the  plate  by  a value  which  corresponds  to  the 
increase  in  the  wing  lift  on  length  dx1>  i.e., 

£ - r.'V- 

Integrating  this  equation  along  the  length  of  the  cord 
0 < x1  < b,  we  obtain  the  wing  lift 


■J 


4Y 


rfx,  « V*  Ml  (*). 


(1) 


where  m(b)  - the  attached  mass  of  the  plate  which  is  determined 
in  the  wing  cross  section  along  the  maximum  span. 


The  value  of  the  transverse  aerodynamic  moment  relative  to  the 
axis  OZ^  which  passes  through  the  apex  of  the  wing  is  calculated 
from  the  formula 


r dy  r 

x (•(*)* -j «(«,)*,}. 


(2) 


The  pressure  drag  (induced  drag)  of  the  wing  with  consideration 
of  the  suction  effect  on  the  leading  edges  is  determined  from  the 
relationship 


Xi-±k«. 


(3) 


Let  us  introduce  into  consideration  instead  of  the  forces  and 
moment  (l)-(3)  the  corresponding  coefficients 

where  S - the  area  of  the  wing  in  a plane ; 
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- the  dynamic  head; 


Then  the  aerodynamic  characteristics  of  a wing  with  small 
aspect  ratio  can  be  presented  in  the  following  form 


If  the  coefficients  of  lift  and  lateral  moment  are  known,  the 
position  of  the  wing's  center  of  pressure  in  relation  to  its  apex 
in  fractions  of  the  root  cord  can  be  calculated  from  the  formula 

(5) 

Formulas  (4)  and  (5)  were  obtained  from  the  general  relation- 
ships of  the  theory  of  a thin  body  and  are  suitable  in  all  cases 
where  the  assumptions  of  this  theory  are  valid.  Let  us  note  two 
special  features  which  follow  from  formulas  (4)  and  (5). 

1.  In  accordance  with  (4)  the  lift  and  moment  coefficients 
are  a linear  function  of  the  angle  of  attack.  The  results  of 
experimental  studies  confirm  the  linear  nature  of  these  relation- 
ships for  wings  of  small  aspect  ratios  to  a < 6°. 

2.  The  lift  coefficient  does  not  depend  on  the  shape  of  the 
wing  in  plane  z1  = i(x1)  and  is  determined  only  by  the  wing  cross 
section  on  the  trailing  edge.  Experimental  data  and  calculations 
in  accordance  with  precise  theories  also  confirm  this  special 
feature  for  wings  with  an  aspect  ratio  X s 1.5. 
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Thus,  despite  limitedness  by  the  case  X 0,  the  theory  of  a 
thin  body  leads  to  qualitatively  correct  results  for  wings  of 
finite  aspect  ratio-  It  can  be  assumed  that  under  certain  con- 
ditions this  theory  will  also  provide  satisfactory  quantitative 
results . 


From  (3)  and  (4)  it  follows  that  the  task  of  determining  the 
aerodynamic  characteristics  of  a wing  of  small  aspect  ratio  in  a 
restricted  flow  is  reduced  to  finding  the  connected  mass  of  a 
plate  close  to  the  surface  of  separation.  To  calculate  the  con- 
nected mass,  we  use  the  method  of  special  features  and  we  distri- 
bute the  vortical  layer  with  continuous  intensity  y(z)  along  the 
length  of  the  plate  within  limits  -l  < z < l.  Obviously,  in  this 
case  y(z)  ■ -y(-z).  We  consider  the  effect  of  the  surface  of 
separation  by  the  mirror  reflection  method  (Fig.  1). 


The  distribution  y(z)  should  satisfy  the  following  condition 
on  the  plate:  a particle  of  the  fluid  which  is  adjacent  to  the 
plate  at  any  point  S should  possess  a constant  vertical  velocity 
V0  (Fig.  1).  Using  this  boundary  condition,  we  obtain  the  integral 
equation  for  the  determination  of  the  unknown  function  y(z) : 


!»*>(*-*,)* 


2*V%. 


(6) 


Here,  the  "plus"  sign  pertains  to  the  movement  of  a wing  beneath 
a free  surface  (hydrofoil),  and  the  "minus"  sign  - to  the  movement 
of  the  wing  close  to  a solid  surface  (ground).  The  values  of  l 
and  h are  determined  in  the  lateral  plane  of  flow  ZOY  depending 
on  coordinate  x1  which  enters  into  equation  (6)  as  a parameter. 

To  calculate  the  lift  and  induced  drag  coefficients  of  the  wing 
the  values  of  l and  h should  be  determined  in  the  cross  section  of 
the  maximum  span. 


z 


We  introduce  the  new  variable  $ from  the  relationship 
l cos  * (with  -Z  c z < Z we  have  ir  < * < 0)  and  we  conduct  the 


replacement  y(z)  = y(l  cos  } ) = y($),  designating  h = h/b . Then 
equation  (6)  Is  written  as  follows: 


± f g«Vt. 

) COT  * — CM  *,  J (COI*  — «•»,)•  + 4A* 


(7) 


Let  us  note  a partial  solution  of  equations  (6)  and  (7)  which 
corresponds  to  the  motion  of  a wing  in  a limitless  fluid  (h  = «>) . 
From  (6)  and  (7)  with  h = 00  we  obtain 

JH??- j-sSSSr*-**. 

v^ience 


7-2V.77ff-?-2V.dg». 


(8) 


For  the  solution  of  equation  (7)  in  the  general  form  with 
h / » we  present  the  unknown  function  y(^)  as  a series 

T ■=  2V*  + JJ  ^ia«in 2n •] . (9) 

in  which  the  first  term  with  AQ  = 1 is  determined  by  the  solution 
of  (8)  and  the  second  term  and  AQ  / 1 characterize  the  additional 
load  which  arises  on  the  plate  from  the  influence  of  the  flow 
boundaries.  Obviously,  with  h •*  «>  coefficient  AQ  -*■  1,  and 
A2n  * 0 (n  - 1,  2 ...). 

In  recording  series  (9)  we  used  the  condition 

% 

If  the  solution  to  (9)  is  known,  then  the  associated  mass  is 
calculated  in  the  following  manner.  The  values  of  the  function 
of  the  velocity  potential  on  the  surface  of  the  plate  with  VQ  ■ 1 
is  determined  with  sin  accuracy  to  a constant  which  subsequently 
is  not  significant: 
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where  the  upper  sign  pertains  to  the  upper  surface  (<f>  ),  and  the 

B 

lower  - to  the  lower  surface  (<J>  ).  The  general  formula  for  the 

H 

determination  of  the  associated  mass  has  the  form  LI] 

« — 

t 

Since,  in  this  case,  on  the  lower  surface  of  the  plate  + 
and  on  the  upper  1,  and,  besides,  dS  = dz, 

s 

i T 

«- -«f{ t. <*><*-  -V  J f- (*>•*«»•**  - 

« 

T 

= 2f /*  j*in»rf§iT(t)sin»rf». 


Substituting  the  expression  y(*)  from  (9)  into  this  formula. 


we  obtain 


» + tP)- 


Thus,  to  determine  the  associated  mass  of  a plate  it  is 
necessary  to  know  the  first  two  coefficients  AQ  and  A 2 of  the 
expansion  (9).  The  remaining  coefficients  A^n  (n  = 2,  3,  4...) 
characterize  the  load  distribution  over  the  span  of  a plate  and 
do  not  have  a direct  influence  on  the  total  aerodynamic  charac- 
teristics of  the  wing. 

With  h * »,  for  the  solution  of  (8)  it  follows  that  AQ  = 1 
and  A2n  * 0;  therefore,  it  accordance  with  (10) 

(11) 

The  aerodynamic  characteristics  of  the  wing  with  consideration 
of  (11)  are  determined  from  formulas  (4)  and  (5) 
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■.  w,  — yUg-« 

1 ~JpW{ 


Here  V - the  aspect  ratio  of  the  wing 
equation  of  the  form  of  a wing  in  a plane. 


Formulas  (12)  are  known  relationships  of  a wing  of  small  aspect 
ratio  in  a limitless  fluid  [2], 


Let  us  move  on  to  the  calculation  of  the  expansion  coefficients 
A2n(  n = 0,  1,  2...).  Equation  (7)  with  consideration  of  (9)  is 
transformed  as  follows: 


where  we  introduce  the  designations 


/(it.  -eo.».y1  + ^y,)- 

T«»*»  £ A*,  (Jtm-l  4-  Jtm+l)  4-  -j  ^ — J to+t). 


After  calculation  of  the  integrals  Jm  and  further  transfor- 
mations of  formula  (13)  we  finally  obtain  the  equations  for  the 
expansion  coefficients  A-  (n  = 0,  1,  2,  ...). 


1.  For  the  motion  of  a wing  close  to  a solid  surface 


2.  For  the  motion  of  a wing  beneath  a free  surface  (hydrofoil) 


+ ••)]  — 1 — jE  4*  [<•» (*,  *0  — eo*2n»J 


The  following  designations  are  introduced  into  (14)  and  (15): 

C0  {h<  ft,)  — — 1 cos  ft,  | — 2 hy\ 

±^fnT^y. 

a-|KK-(sin,»,  + 4 ft; 

5-4.Kv,+  *in»ft.  + 4A*; 

V = Y (sin*ft,  + 41)1  + 16  A*  cos* 

t = == ; f 2,  (/)  = cos  (2n  arc  cos /), 

M'  + ft* 

where  T2  (t)  " the  Chebyshev  polynomial  type  I. 

If  in  relationships  (14)  or  (15)  we  record  a number  of  values 

of  } , we  obtain  systems  of  algebraic  equations  for  the  determin- 
s 

ation  of  expansion  coefficients  A2n(n  = 0,  1 ...). 

Here,  the  number  of  recorded  points  $ . (i  = 1,  2 ...)  determines 
the  number  of  equations  in  systems  (14)  and  (15)  and,  consequently, 
the  number  of  unknown  coefficients  A2n  (n  = 0,  1,  2,  ...,  i - 1). 

The  remaining  coefficients  (n  = i,  i + 1,  ...)  should  be  taken  as 
equal  to  zero. 

In  accordance  with  (14)  and  (15),  calculations  of  the  coef- 
ficients were  conducted  with  the  fixing  of  1 point  (approximation 
I),  three  (approximation  II),  and  five  (approximation  III). 

Some  of  the  calculations  for  the  motion  of  a wing  close  to  a 
solid  surface  are  presented  in  Fig.  2.  From  the  known  coefficients 
Aq  and  A2  and  formula  (10),  we  can  find  the  associated  mass  of  the 
plate,  and  from  formulas  (4)  and  (5)  - the  aerodynamic  character- 
istics of  the  wing. 

Figure  3 presents  the  values  of  the  derivative  of  the  lift 
coefficient  in  accordance  with  angle  of  attack  CJ,  referred  to  the 
value  C^m  in  a limitless  flow  (12).  Here,  the  dots  designate 
the  mean  values  of  CJ,  which  were  obtained  in  accordance  with  the 
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theory  of  a vortical  lifting  surface  for  rectangular  wings  with 
A = 0.4-1. 6 [4]. 

An  analysis  of  the  calculations  which  were  conducted  shows 
that  satisfactory  results  right  up  to  h 2.  0.1  can  be  obtained 
already  in  approximation  III  for  a wing  close  to  a solid  surface 
and,  in  II  - for  a hydrofoil. 

If  we  restrict  ourselves  to  values  of  h > 0.4  (wing  close  to 
the  ground)  and  h > 0.25  (hydrofoil),  the  aerodynamic  character- 
istics of  a wing  with  a small  aspect  angle  can  be  determined  from 
approximation  I.  In  this  case  simple  analytical  relations  occur 

which  follow  from  (14)  and  (15)  with  $ = ir/2  and  A„  = 0 (n  = 1, 

s cL  n 

2,  3 ...). 


1.  For  a wing  close  to  a solid  surface 


Since 


Ki-nfr 

£ * 


m <=  *p/* 


VT+fr 


Then 


C'  3 


2.  For  a hydrofoil 

= - 


: — sr" 

a-77T3T 


m = *p  /*• 


j— 


Vi+< 


1— 


(17) 


The  calculations  from  formulas  (16)  and  (17)  are  shown  in 
Fig.  2 by  the  solid  lines. 
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we  replace  2 by 


Let  us  note  that  If  in  the  expression 
a larger  number  (for  example  3)»  formulas  (16)  and  (17)  provide 
a satisfactory  agreement  with  calculations  in  accordance  with 
approximation  III  right  up  to  h i 0.10. 


Thus,  we  can  recommend  the  following  approximate  formulas  for 
the  calculation  of  the  associated  mass  of  a plate  in  a restricted 
flow  with  h > 0.10. 


1.  Close  to  a solid  surface 

(18) 

2.  Under  a free  surface 

""7  "'rf  (19) 

Formulas  (18)  and  (19)  are  convenient  with  a calculation  of 
moment  characteristics  and  the  position  of  the  wings'  center  of 
pressure.  For  example,  for  a wing  close  to  a solid  surface  from 
(5)  and  (18)  we  find 


jru,>yz 


CTOLMD 


*u») 


dMt. 


Before  the  integral  h = h/Z(b),  h - the  position  of  the  trailing 
edge  of  the  wing  in  relation  to  the  interfaces,  Z(b)  - the  wing 
semispan. 


In  the  integrand  the  values  h and  l depend  on  the  coordinates 
x^.  With  small  angles  of  attack 

*<*«>  /(»)  [1+ rJjjd -*,>•]. 

In  the  assumptions  of  the  theory  of  a thin  body  , under 

certain  conditions  the  second  term  in  the  right  side  can  also  have 
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the  order  of  the  first  term;  therefore,  the  position  of  the  center 
of  pressure  on  a flat  wing  close  to  a surface  of  separation  depends 
on  the  angles  of  attack. 


Let  us  limit  ourselves  to  the  case  a -*•  0;  then  - Z(b)h 
= const  and  the  position  of  the  center  of  pressure  is  determined 
from  the  formula 

**- 1 - f * V 


(20) 


where  T(x,)  —7^  - the  relative  local  wing  semispan; 


Let  us  examine  a family  of  wings  in  which  the  form  in  a plane 
is  given  by  the  equation  i(x^)  = x-™. 


The  exponent  m varies  within  limits  of  0 < m < ®,  With  m = 1 
we  obtain  a delta  wing  and  with  m 0 - rectangular. 


In  accordance  with  formula  (20)  we  will  have 

1 


*ry*r+1 


(21) 


If  we  accomplish  such  conversions  for  a hydrofoil,  then  from 
(5)  ana  (19)  with  a -*■  0 it  follows  that 


■ -(»-7rbp)  jiWXf+V-* 


(22) 


For  a family  of  wings  T(x^)  = x™ 

, L * \frv*rr  r 


(23) 


11 


With  h 


S 00 


(the  flow  of  a limitless  fluid)  from  (21)  and  (23) 


t 


I 


(24) 


This  result  can  also  be  determined  from  formula  (12). 


In  the  general  case  the  integrals  in  formulas  (21)  and  (23) 
are  expressed  by  hypergeometric  functions  [3].  Figure  4 presents 
the  calculations  of  for  wings  with  m = 1/2,  l,i(solId 

lines  - for  a wing  close  to  a solid  surface,  dashed  lines  - for 
a hydrofoil). 

An  analysis  of  the  results  which  have  been  obtained  shows  that 
with  the  approach  of  a wing  to  a solid  surface  (ground)  the  center 
of  pressure  is  shifted  toward  the  trailing  edge,  in  which  regard 
this  displacement  is  greater  the  smaller  the  exponent  m.  In 
particular,  the  greatest  displacement  should  be  observed  for  a 
rectangular  wing  (m  0). 

The  effect  of  the  boundary  of  a free  surface  for  a hydrofoil 
Is  opposite  to  the  effect  of  the  ground. 
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Key:  (1)  Approximation 

(2)  Solid  surface. 


Key:  (1)  Approximation;  (2) 

According  to;  (3)  Solid 
surface;  (4)  Free  surface. 
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